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Abstract. Consider a compact manifold M with boundary 8M endowed with 
a Riemannian metric g and a magnetic field Q. Given a point and direction 
of entry at the boundary, the scattering relation E determines the point and 
direction of exit of a particle of unit charge, mass, and energy. In this paper we 
show that a magnetic system (M, 8M, g, f2) that is known to be real-analytic 
and that satisfies some mild restrictions on conjugate points is uniquely deter- 
mined up to a natural equivalence by E. In the case that the magnetic field 
Q is taken to be zero, this gives a new rigidity result in Riemannian geometry 
that is more general than related results in the literature. 



1. An introduction including the results PROVED 

1.1. Magnetic systems and scattering data. Let M be a compact, 71-dimensional, 
Riemannian manifold with boundary dM and metric g. Let tt : T*M — > M be the 
natural projection of the cotangent bundle to M given by (x,^) >—>■ x. The Hamil- 
tonian flow of the function H(x, £) = with respect to the standard symplectic 
structure of T*M is the geodesic flow. Given a closed 2-form Q on M, there is a 
corresponding magnetic geodesic flow which is determined by the same Hamiltonian 
function, but through the modified symplectic 2-form 

uj = dx A dt; + 7T*0. 

The integral curves of this flow, projected to M by ir, are the magnetic geodesies. 
The level sets of H are preserved along the flow and that implies that the magnetic 
geodesies have constant speed. In this paper, we shall restrict our attention to 
the energy level H = 1/2, the unit co-sphere bundle S*M. The isomorphism 
T*M — > TM induced by the metric g transforms the magnetic flow to a flow on 
SM, the sphere bundle in TM. Let X and denote the vector fields on SM 
generating the standard geodesic flow and magnetic geodesic flow, respectively. 
They are related by the equation 

Xpfav) = X(x, v) + Y{x)v, 

where Y is an antisymmetric (1,1) tensor on M satisfying 

(Y(x)v,w) g = Cl(v,w). 

Let $' : SM — > SM denote the magnetic geodesic flow. Given a pair (x, 8) € SM, 
let 7^.9 or j x .e(t) denote the magnetic geodesic with initial point and vector (x,9). 
We also use the magnetic exponential map 

cxp^ : TM -> M, 

defined by exp^(t#) = jx,e{t) for unit length 6. This map is C 1 everywhere and 
smooth away from 9 = 0. 

1 



2 



P. HERREROS AND J. VARGO 



Definition 1. A 4-tuple (M,dM,g,il) is called a magnetic system. Two such 
systems (Mi,dMi,gi,Cli), i = 1,2 with the same boundary dM\ = dMi are said to 
be equivalent if there exists a diffeomorphism (p : Mi — > Mi such that: 

<P*92 = 0i, <P*^2 = Hi, <y9| aA f = id. 
The diffeomorphism ip is called a magnetic equivalence. 

In this paper, the term analytic shall mean real-analytic. We define a magnetic 
system (M, DM, y, fi) to be analytic if there exists a magnetic system (M', DM', g', CI') 
with M' a real-analytic open manifold containing M such that: 

(i) g'\ M = 9 and fl'\ M = 0,, 

(ii) g' and fl' are real-analytic on M', 

(iii) 9M is a real- analytic submanifold of M'. 

The extensions (g' , 17') are uniquely determined from (y, f2) by analytic continuation 
and shall henceforth be denoted without the prime. That is, (y, fi) shall denote its 
own extension to M' . 

The scattering magnetic rigidity question is whether the equivalence class of a 
magnetic system is determined by the scattering of its magnetic geodesies. Given 
a point of entry x £ dM and an initial direction 8 £ S X M, the scattering data tells 
all future points of contact j X) g n dM as well as the directions of motion at those 
points. More precisely: 

Definition 2. A pair of vectors (v, w) £ dSM x dSM belongs to the scattering 
relation E if for some t >Q, to = 

After properly defining all the terms and introducing some minor restrictions on 
conjugate points, it is shown that two analytic magnetic systems are equivalent if 
they have the same scattering data. 

This problem descends from the boundary rigidity problem in which one seeks 
to recover a metric g up to isometry from the boundary distance function p g : 
dM x dM -> M, defined by: 

p g {x,y) = inf{length(7) : 7 C Mis a curve joining x toy}. 

A particular Riemannian manifold with boundary (M, dM, g) is called boundary 
rigid if all other such manifolds with the same boundary distance function are in 
the same isometry class. The known examples of manifolds that are not boundary 
rigid cither have trapped geodesies (geodesies that never leave the manifold) or con- 
jugate points. The most general class of manifolds conjectured to be boundary rigid 
are those with the SGM property. This property is detectable from the boundary 
and is equivalent to the two conditions that there are no trapped geodesies and that 
all geodesic segments are strong length minimizers [CI j . SGM manifolds homeo- 
morphic to a ball with a strictly convex boundary are called simple. It has been 
conjectured that all simple manifolds are boundary rigid (e.g. |Mij). but so far, 
such a general statement has only been confirmed in dimension 2 in work by Pestov 
and Uhlmann PUJ. In higher dimensions, Gromov showed that subregions of Eu- 
clidean space are boundary rigid [Gr| . Other spaces known to be boundary rigid 
are subregions of boundary rigid spaces, subregions of a hemisphere [Mi] , simple 
symmetric spaces of negative curvature }BCG| , and SGM subregions of spaces with 
a parallel vector field [CKj . Recently, it was shown by Burago and Ivanov that 
metrics sufficiently close in the C°° topology to the Euclidean metric on a convex 



SCATTERING RIGIDITY FOR ANALYTIC MANIFOLDS WITH A MAGNETIC FIELD 3 

domain are boundary rigid [BIj . The last result implies, for example, that for any 
point on any Riemannian manifold, a sufficiently small neighborhood of that point 
is boundary rigid. 

In applications, this problem is called travel time tomography and goes back at 
least as far as the early 20th century, when it was proposed to measure seismic 
waves to reconstruct an image of the inner structure of the earth. Besides being 
used widely by geophysicists, travel time tomography is used by occanographers and 
atmospheric scientists, and has potential applications in medical imaging. In some 
applications, it is natural to assume that the metric is isotropic; that is, conformal 
to the Euclidean. For metrics in a given conformal class, uniqueness was first shown 
for simple metrics by Mukhometov [Muj and later by Croke for SGM metrics jClj . 

The SGM and simplicity conditions are restrictive geometric conditions which 
are not typically satisfied in applications. In more general circumstances, one can 
instead consider the lens data, which consists of the scattering data together with 
the travel times of the geodesies. It was shown by Michel that the lens data is 
equivalent to the boundary distance function for simple manifolds |Mij . If trapped 
geodesies are allowed, then counterexamples to lens rigidity exist |CKj . Croke 
showed that if a manifold is lens rigid, then so is a finite quotient of that manifold 
[C2j . In [V] lens rigidity was shown if the metrics are a priori assumed to be 
real analytic. Stefanov and Uhlmann proved lens rigidity for metrics that are a 
priori known to be close to a given base metric taken from a set of generic metrics 
including the real analytic ones |SU2j . 

The magnetic rigidity problem is related to the question of whether there are 
more general families of curves whose defining parameters can be recovered from 
boundary measurements. Magnetic geodesies and the magnetic flow were first con- 
sidered by V.I. Arnold [Arlj and D.V. Anosov and Y.G. Sinai [AS] , Since then 
it has been studied from several approaches (e.g. |CMP( IGi( |Grog| INT| IPPj ) and 
their boundary rigidity, in particular, has been studied in [DPSU] IHej . 

One motivation for magnetic boundary rigidity is that it has potential appli- 
cations to problems of geometry. In the case of surfaces, magnetic geodesies are 
related to the study of curves of constant geodesic curvature. This relation has been 
used to study the existence of closed curves with prescribed geodesic curvature (see 
e.g |Ar2( ILe[ ISchj ) . Magnetic boundary rigidity was used in [Hej to show that a 
surface of constant curvature cannot be modified in a small region while keeping 
all the curves of a fixed constant geodesic curvature closed. 

Let l : dM M be inclusion. Since a magnetic equivalence fixes the points 
of dM, the pullbacks i* g and t*fl are preserved. However the scattering relation 
£ is not preserved because its definition depends on the set dSM which is defined 
through g\dM- For this reason we reparametrize £ through the natural diffcomor- 
phism 

A : dSM -> dSM = {{w, r) G TdM x K : \w\ 2 g + r 2 = 1}. 

A is given by v >— > (v T ,(v,v) g ), where v T denotes the orthogonal projection of 
v G dSM to TdM, and v denotes the inward-pointing unit normal vector. The 
reparametrized scattering relation, £ C dSM x dSM, is defined to be £ = (A x 
A)(£). Given an equivalence class of magnetic systems {(M, dM, g, fi)}, we shall 
call (dM, i*g, £) the associated scattering data. 
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An equivalent way of parametrizing the scattering data is given through boundary 
normal coordinates. Given xo £ dM , there exists a neighborhood N C dM of xq 
and a number e > such that the map 

exp,, : N x (-e, e) -> M' 

given by exp 1/ (x / , x n ) = exp x ,(x n p) is an injective immersion. Here exp denotes 
the exponential map corresponding to the Ricmannian metric g, not the magnetic 
exponential exp M . Moreover, if x a : 1 < a < n , — 1 are local coordinates for dM, 
then the metric tensor has the form 

g = g a p(x)dx a dx f3 + {dx n f . 

(x 1 ', x n ), as local coordinates for M near xq are called boundary normal coordinates, 
or semigeodesic coordinates. Note that the curves given by x' = const are geodesies 
(not magnetic geodesies) of (M, g) normal to dM. 

Since dM is compact, it is possible to extend the boundary normal coordinates 
to a collared neighborhood of the entire boundary. Indeed, if e is sufficiently small, 
then we obtain a diffeomorphism 

(1.1) exp„ : dM x (-e, e) V, 

where V C M' consists of the points whose distance to dM is less than e. For a 
detailed proof, see [V] . It shall be convenient to set 

M 1 = M U V. 

Finally, when discussing two magnetic systems with the same boundary, we shall 
assume that the value of e is chosen sufficiently small so that (jl.lj) is a valid diffeo- 
morphism for both. 

Note that the mapping A is equal to the pullback of vectors at the boundary by 
cxp^. 

A = exp* \qsm- 

Definition 3. (M,dM,g,Sl) is called non-trapping if for all (x,0) £ SM, there 
exist T_ < and T + > such that r ) x ,e{T±) G M' \ M; that is, every magnetic 
geodesic must eventually leave M in both directions. 

In this paper, we shall assume that all magnetic systems are non-trapping. Let 
v = (x, 9) G dSM, and let t{v) = T be defined by 

T = mf{t>0: lv (t) <£ M}. 

Then j v (T) and j' v (T) are respectively called the terminal point and direction of 
7„ in M . The mapping 

t : dSM -» R 

is called the travel time map of M. The corresponding map £ : dSM — > M given by 
£ = £ o A -1 is called the travel time data of the equivalence class of (M, dM, g, f2). 

Definition 4. We say that v satisfies condition A if there are no points in dM n 
{jv{t) ■ < t < T} that are conjugate to x along j v . In this paper, conjugate 
shall always mean with respect to the magnetic flow. We shall say that an analytic 
magnetic system (M,dM, g,fl) satisfies condition A if there exists at least one 
vector in each connected component of dSM that satisfies condition A. 
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Let Q be the collection of magnetic geodesic segments in M whose endpoints lie 
in dM and do not intersect dM at any other point . Let 7 G Q be an element with 
initial vector v G SM, travel time T, and final point y — j(T). 

Definition 5. We say that 7 satisfies condition B if there are no points on 7 that 
are conjugate of order n — 1 to y. We shall say that an analytic magnetic system 
(Af, dM, g,fl) satisfies condition B if the set of magnetic geodesic segments in Q 
satisfying condition B is dense in Q . 

Recall that if x and y belong to M and y = exp^(t>), then x and y are conjugate 
if d v exp%(v) is singular. The order, or multiplicity, of conjugacy is equal to the 
dimension of the kernel with n — 1 being the largest possible. For example, two 
antipodal points on S n (with = 0) are conjugate of order n — 1. 

Lemma 1. Condition A is open in dSM. Condition B is open in Q . 

Proof. The openness of condition B is clear. Let vo G dSM satisfy condition A. 
Let K = j Vo n dM. By compactness, for any open neighborhood N of K, there 
exists an open U about vo in SM' such that for all v G U, 7„ n dM C N. 

Since xq, the basepoint of vq, is not conjugate to any point of K, we may choose 
N and the corresponding neighborhood U sufficiently small so that for each v G U, 
its basepoint x will not be conjugate to any point of j v P\N. But since j v (ldM C N, 
we find that x will not be conjugate to any point of j v n dM. □ 

Theorem 1. Let [Mi,dM,gi,0,{), i = l,2 be two non-trapping, analytic magnetic 
systems satisfying condition B and assume that (Mi, dM, gi,Qi) satisfies condition 
A. Then if the two systems have the same scattering data, they must be equivalent. 

The first step of the proof is to show that the mapping ipo : V\ — > V2 given by 

(1-2) <Po = exp„ 2 o exp^ 1 

is a magnetic equivalence fixing the points of dM. To do this, we must use condition 
A to show that the coefficients of (g, f2) are uniquely determined by the scattering 
data in boundary normal coordinates. The main step is to prove the following 
theorem. 

Theorem 2. Let (Mi,dM, 5i, Oj), i = 1,2 be non-trapping, analytic magnetic sys- 
tems with the same scattering data, and assume that vo = (xo,@o) G SdM satisfies 
condition A in Mi. Then at the point xq, the coefficients of (gi, Oi) and (52,^2), 
expressed in boundary normal coordinates, have the same jets at Xq. That is, their 
values at xq are the same, and the values of all their derivatives of all orders at xq 
are the same. 

This theorem is a direct generalization of a result by Plamen Stefanov and Gun- 
thcr Uhlmann for Ricmannian metrics (with no magnetic fields). We extend their 
result by allowing magnetic fields and by only assuming the scattering data and 
no magnetic analog of length data. In order to generalize the proof, some extra 
condition is needed for the case when the magnetic geodesies near 7„ are short. 
For this paper the convenient extra assumption is to assume that the system be 
analytic. However it should be noted that this is not at all an essential condition 
and could be easily be replaced by a certain weak convexity condition. See the 
remark after the proof for details. 
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A similar theorem was proved for magnetic systems by |DPSU| in which they 
assume strict convexity of the boundary with respect to magnetic geodesies and 
take a certain action function as the given data. We use a similar action function 
to prove this theorem, but do not otherwise use the same techniques since we allow 
for the possibility that j Vo be a long magnetic geodesic. Rather we follow Stefanov 
and Uhlmann's proof, the idea of which is to consider a function p[x, y) which gives 
travel times between pairs of points. Fixing y, the function satisfies the eikonal 
equation: 

g ll (x)(d XtP )(d X] p) = 1. 

Using this equation, the scattering data tells us just enough about the derivatives 
of p to reconstruct the jet of <?. The proof here follows the same program, but with 
modifications to account for magnetic fields and the lack of explicit travel time data. 
We construct an analogous action function p(x, y) which satisfies an eikonal-type 
equation (equation (|2.4|) ) whose coefficients come from g and O. We then use the 
scattering data to infer enough information about the derivatives of p to recover 
the jets of g and f2. 

The next step in the proof of Theorem Q] is the proof of the following theorem. 

Theorem 3. Let (Mi,dM,gi,£li), i = l,2 be two non-trapping, analytic magnetic 
systems that satisfy condition B and have the same scattering data. If tfo : V\ — > V2 
is an equivalence, then corresponding elements of Q\ and Q2 have the same travel 
times. In particular, the two systems have the same travel time data I. 

Elements of Q\ and Qi are corresponding if their initial and final vectors are 
related by ipo. The next theorem finishes the proof of [TJ 

Theorem 4. Let (Mi, DM, g^fij), i = 1,2 be two non-trapping, analytic magnetic 
systems with the same scattering and travel time data. If ipo : V\ — > V2 is an 
equivalence, then ifo extends to an equivalence tp : M[ — > M^. 

2. Theorem [2] Recovery of the magnetic system in a neighborhood of 

the boundary 

2.1. Defining an action function on magnetic geodesic segments. The first 
step shall be to construct an action function which serves the role of the travel time 
function used by Stefanov and Uhlmann. Subscripts shall be omitted as long as all 
statements apply equally to both manifolds. We consider our magnetic system to 
be extended to an open manifold M 1 containing M. Let 7 = {7(f) : < t < T} be a 
magnetic geodesic segment in M' . For our proof, we shall need a 1-form £, defined 
in a neighborhood of 7, satisfying d( = Q. By the Poincare lemma, such a 1-form 
exists if 7 has no self-intersections. If 7 does intersect itself, then we circumvent the 
topological problem through the construction described in the following paragraphs. 
Let 7 denote the segment in M. n given by 

7(t) = (0, . . . , 0, t) for < t < T. 

Then there exists a neighborhood U C M. n of 7 and an immersion tp : U — > M' that 
satisfies ip(0,t) = ipo^f(t) = 7(f). For example, we can construct such a map in the 
following way. Extend 7 to an open interval containing [0, T] and choose vectors 
Xi(t) along 7(i) such that (X\(t), . . . , X„_i(t), 7'(f)) arc linearly independent for 
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all t. Then define 

i>{x\...,x n -\x n )=vq> l(xn) (^x'X^x^ . 

The topology of U is trivial so by the Poincarc lemma, there exists £ such that 
dQ = tp*n. We call £ a magnetic potential in a neighborhood of 7. 

Let i: be coordinates for U , and (i, £) the corresponding natural coordinates for 
T*U . Using ip we can pull back structures from M' to U. In particular, wc obtain 
an exact symplcctic form 

Q = dx A d£ + ip*Cl = -d (idx - C) 

and a Hamiltonian function H = Together, these generate Hamiltonian 

curves whose projections to U are the magnetic geodesies of the magnetic system 

{u,rg,pn). 

Wc parametrize the magnetic geodesic segments in M' near 7 by their initial 
vector v £ SM' and their length r. 

(v, t) h-> 7„ jT = {7t,(t) : < i < t}. 

Let (i>o,7o) correspond to the original segment 7 about which £ was constructed. 
For (v,t) in some neighborhood M. about (vq,tq), the magnetic geodesic segments 
7„. T in M' can be uniquely pulled back via ip to magnetic geodesic segments >y VT 
in U. Let c„ iT denote the corresponding Hamiltonian curve in T*U. 
Letting J\f denote the neighborhood of magnetic geodesic segments 

N = {lv,r ■ (v,t) e M}, 
we define the action functional : J\f — > K by 

A c [7v,r] = r- C- 

In the case that fi is exact, £ can be defined on all of M' so that is well-defined 
for all magnetic geodesies. But in general this is not possible so the action can 
only be defined near one fixed curve. In cither case, A^ depends on the choice of £ 
which is not uniquely determined by the magnetic system. 

Proposition 1. 

A C [ 7v> r]=/ (£dx-7T*0 

Proof. Along the curve c„ )T , we have ^dx = |£|gdt = dt. Therefore, noting that 
Jv (t) = 7T c v (t) , the integral on the right is equal to 

T-f n*C = r- [ C 



□ 
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2.2. The first variation of Af. Let J s (t) : — e < s < e, < t < t s be a smooth 
1-parameter family of unit-speed magnetic geodesies in TV. Correspondingly, there 
is a smooth 1-paramctcr family of curves c s (t) in T*U , and by the lemma, we have 



tp(a) = A c [ 7s ] = / (fdx-7r*C). 



Let a(s) be the curve c s (0) and 6(s) be the curve c s (r s ). Now let < \h\ < e, and 
consider the surface a with paramctrization (s,t) i— > c s (i) : < s < ft, < t < t s . 
According to Stokes' theorem, 

' (idx- 7T*C) = [(diAdx- 7T*ft). 



CO 



It follows that 



The surface integral on the right-hand side is equal to 0. To show this we note 
that the 2-form being integrated is — u), the symplectic form. We also note that 
the Hamiltonian vector field = JjC s (t) lies tangent to a . By definition, 

w(X & , -) = dH. 

Since a is contained in the level surface H = 1/2, we conclude that the pull back 
of l> to a must be 0. Hence we obtain 
<f{h) - y(0) _ 1 
h h 
Taking the limit as h — > 0, we find 



(£di-7T*0. 



^(0) = (^- 7 r*C)| 6 '/(°o ) ) 

c s (t) = (7s(*),^ s (i)) where ^ s (i) is the covector corresponding to 9t7 s (<) by g. 
Therefore 

(2.i) ^(o) = [(5 t 7o, -) 9 -c]\:: b X\ 

In terms of the family of magnetic geodesies {■j s (t)} s , 7r*a'(0) and 7r*b'(0) are just 
the variations in the initial and terminal point, respectively, of the curves. 

Now suppose that 7 = 7„ 0jTo h as endpoints xq and yo which are not conjugate 
to each other. Then for (x,y) in a sufficiently small neighborhood around (xo,j/o)j 
there exists a magnetic geodesic "f XtV joining x to y which is smoothly dependent on 
x and y. Let 6 = 6 x y G T X M' be its initial vector and define the smooth function 

p(x,y) = A c [7 x ,y]. 

Let w £ T X M' . If we fix y and take the differential of p with respect to x, then by 
(12.11). we obtain: 



(2.2) (d x p,w) = -(6,w} g + ({,w) 
Solving for 8, we find 

(2.3) -(9,-} g =d x p-C 

In particular since 7 is a unit-speed curve, 9 has unit length, so 

\d x p-(\ 2 q = l. 
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Written another way, 

(2.4) 9 ij (x)(d xi p - Q(x))(d x3 p- Cj(x)) = 1. 

This is the eikonal-type equation that we need. 

2.3. Choosing compatible magnetic potentials. Now we consider again our 
two manifolds M\ and Mi with their respective metrics and magnetic fields and 
extensions M[ and M' 2 . We let v a = (x ,0 ) <G SdM be a vector satisfying the 
hypothesis of the theorem. For each i = 1, 2, let ji C M[ be the magnetic geodesic 
with initial vector vq, and let Q be an arbitrarily chosen magnetic potential in a 
neighborhood of 7,. The next lemma shows that the two potentials can be made 
compatible with each other in a certain sense. Recall that (x' , x n ) denote boundary 
normal coordinates near xq. 

Lemma 2. // the two magnetic systems have the same scattering data, then there 
exist magnetic potentials (1 , C2 defined in neighborhoods of 71 and 72 respectively 
which satisfy 

hCi = h&i and 

(CiA-> = (Ca,a ai »> ) 

in a neighborhood of the point Xq . 

Proof. Starting with arbitrary ^ and £2, note that 

^(<-iCi _ hb) = '1^1 - = 0. 

Therefore, there exists a smooth function fo defined in a neighborhood of xo in dM 
such that 

t>l& - '-2C2 = df . 

We extend fo in a neighborhood of dM according to the first order PDE 

d x nf=(( 1 ,d x n)-(( 2 ,d xn ), 
f\dM = fo- 

Away from the boundary we apply a cutoff function so that / is well-defined over a 
neighborhood of 71. Letting ([ = (1 +df, wc find that and (2 satisfy the requisite 
conditions. □ 

2.4. Proof of Theorem [2j For the remainder of the proof, subscripts distin- 
guishing the two manifolds shall be omitted. We are given a magnetic system 

(M, dM, g, f2) with known scattering data (dM, i*g, t*fJ, S) and extension (M' , dM, g, O). 
We will be working in boundary normal coordinates near a point xo £ dM. In these 
coordinates, we may regard g\dM as known since L*g is included in the scattering 
data. 

We start with a lemma. Let vq € dSM be a unit vector at the boundary, and let 
C be a magnetic potential in a neighborhood of j Vo . Let yo be a point on 7„ fl dM 
with yo 7^ xo and assume that xq and yo are not conjugate along j VQ . Then there 
exist respective neighborhoods V, W about xo, yo in the extended manifold M' 
such that for any x G V, y &W, there is a magnetic geodesic 7^ connecting them. 
It depends smoothly on the endpoints and is equal to j VQ \[ Xo ,y ] if (x,y) = (xo,yo)- 
As discussed above, the function p(x, y) = A^[7 xy ] is smooth with derivative given 
in equation (|2.2p . 
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According to lemma[2j we may regard i*£ as known data. What's more, we may 
regard f (d x n ) as known in a neighborhood of xo even off the boundary. In particular, 
i*C and C(dx")\dM together give us C|flAf m boundary normal coordinates. If we 
let £ denote the vector field corresponding to £ via g, then that is also known on 
dM. What we need for the theorem is the full jet of g and C at dM . But since we 
already know g and C on dM , what remains is to recover their derivatives in the 
normal direction. That is, we need: 

9^g\ dM , d% n (\dM, A: = 1,2,... 

Actually, what we shall get are these normal derivatives on a half-neighborhood of 
.To in dM. This is enough to apply tangential derivatives to get the full jets at xq. 

Lemma 3. Let v G dSM be a unit vector over the point x G dM which points 
strictly inwards, i.e. (v ,i>)g > 0. Let y be the first point of intersection between 
7„ and dM . Then if xq and go are not conjugate, the scattering data uniquely 
determine d x p(x,yo) for x G dM sufficiently close to xq. 

Proof. The existence of go is guaranteed by the fact that Vq is inward pointing 
and that our systems are non-trapping. Let r > denote the travel time of the 
magnetic geodesic from xq to yo- Then 7„ Q (r) = yo. Let wq = Y v (t), the vector 
tangent to the curve at yo. It should be noted that both yo and wq are determined 
by the scattering relation. For small e > 0, let 

C c = {w G S yo M : (w, w ) g > 1 - e, (w, v) g < 0}. 

Since go is the first point of intersection of 7„ with dM, it is clear that {wq, v) g < 
with equality if an only if wq is tangent to dM. Therefore, C e is either an open 
cone or a half-open cone about wq, depending on whether wo is tangent to dM. 
Note that in the case that it is only half-open, wq is a limit point, but not actually 
a member of C € . 

The elements of C e are strictly outward pointing, so they are terminal vectors of 
magnetic geodesies. Since vo is transverse to dM, we see that {(v, w) G £ : w G C e } 
gives a one-to-one correspondence between vectors v near vq with vectors w G C e . 
Let k denote the composition of the corresponding maps 

C £ -> dSM -» dM, 

where n(v) is projection to the base point. Since xq and yo are not conjugate, 
k is a diffeomorphism. Also, it is determined by the scattering data. Its image, 
which we will denote by H, is either a half- neighbor hood or full neighborhood of 
xq depending on whether C e is a half-cone or full cone about wq. 

The function p(x, yo) is well-defined for x G M' in an open neighborhood about 
xq. We may assume that Ti. is taken small enough to lie inside this neighborhood. 
Then from the previous discussion, we see that for x G H G dM, there exists a 
unique vector v x G S X M, close to ^o, such that y lies on "f Vx . Indeed v x is a point 
in the image of the map C e — > dSM. By equation (|2.2| . 

V x p(x,y ) =( x -v x . 

The terms on the right are uniquely determined, since v x is known from the scatter- 
ing data, and C is known on dM. Therefore, since g\dM is known, we find d x p[x, yo) 
at points x G TL by lowering the index of \7 x p(x, yo). □ 
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Proof of Theorem^ Let vo G SdM satisfy condition A. By lemma HJ the vector 

(2.5) v s (xq) = cos(s)wo + sin(s)i/ 

will also satisfy condition A for sufRently small s > 0. We let y s denote the first 
point of intersection between J v j Xo ) an d dM, and we let w s be the tangent to the 
magnetic geodesic at y s . By compactness, there is a sequence s m — > such that 
(y Sml w Sm ) converges to a vector (y*,u>*) tangent to 7„ and lying over a point of 
dM. We split the proof of the theorem into the two possible cases that u>* = vo or 
to* =/= Vo which correspond to whether the length of the magnetic geodesic j Vo from 
xq to y* has length or not. 
Case I: w* ^ vq. 

Let p m (x) denote the function p[x,y Sm ). Let TL m be the corresponding (half)- 
ncighborhood of xq referred to in the previous lemma. For x € M' nearby xq, let 
v Sm (x) be the initial vector of the magnetic geodesic connecting x to y Sm , and let 
v Sm be the corresponding covector. 

Vs m {x) = ( x - d x p m {x). 

For x £ H-m, both v Sm (x) and v Sm (x) are uniquely determined by the scattering 
data. 

We make the following conventions of notation. Near xq our boundary normal 
coordinates shall be 

{(x a ,x n ) : a=l,...,n-l}. 

A subscript or superscript of a or n shall mean the corresponding tensor component: 
e.g. d n = dx™, Ca = C(da), etc. By equation (|2.4[) . we have in boundary normal 
coordinates: 

g al3 ((a - d a p m )(( - dfiPm) + (On - d n p m ) 2 = 1. 

or 

This equation is valid for all x which are sufficiently close to xq in the extended 
manifold M' . For x £ Tt m , the coefficients of the equation are known. We shall 
find the normal derivatives of g and £ by successively applying d n to this equation 
and solving for all unknown terms at each step. 

For the first step, we apply d n and obtain the equation: 

(2.6) (d n g al3 )(( a - d a p m )(C P - d p m ) + 2g aP (d n ( a - 9„9 Q Pm)(C/3 - dppm) 

+2((„ - d n p m )(d n ( n - d 2 p m ) = 0. 
Splitting the second term and rearranging yields 

(2.7) {d n g ap ){v Sm ) a (v Sm )p + 2(g a ^d n C a )(v Sm ) p = 

2g al3 (v Sm )pd n d a p m - 2(v Sm ) n (dn(n - dlp m )- 
Next, we let to — > oo. On the left side we obtain: 

(d n g a P)(v ) a (vo)p + 2(3 Q/3 d n Co)(«o)/3. 

For fixed x = xq and variable v E S* Q dM, this expression can be regarded as a 
non-homogeneous quadratic function F(v). We wish to show that at x = xq, the 
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value F(vq) is uniquely determined. The fact that condition A is open would then 
imply that the values F(v) are uniquely determined for all unit-length v in an open 
cone. But that implies that the coefficients of F are uniquely determined. Hence 
d n g a (xo) and d n £ a (xo) would be uniquely determined. By again appealing to the 
fact that condition A is open in dSM, we could conclude that d n g a P and d n C,a are 
in fact uniquely determined in a neighborhood of xq in DM. 

To show that F(vq) is uniquely determined, we consider the right side of (|2.7j) 
at x = Xq as 771 — > oo. In the second term, d n Q n is independent of m and 

(2-8) dlp m (x ) = dlp(x ,y Sm ) -> d*p(x ,y*). 

Since (v Sm ) n = sin(s m ) — ► 0, we conclude that the limit of the second term is 0: 

(2-9) 2(v Sm ) n (d n ( n - d 2 p s J ^ 0. 

It only remains to show that the limit of the first term of the righthand side is 
uniquely determined at x$. For each m, d n p m (x) is known in TL m . Hence d n d apm = 
d a d n Pm is also known in TL m and in particular at xq and in the limit m oo. Thus, 
d n g a/3 and d n (, a are uniquely determined for x £ DM near xq. 

Before calculating higher derivatives, we return to equation (|2.7p . All coefficients 
except those in the far right hand term are known in x € Ti. m - Therefore, by solving 
the equation, we find that d n (, n —d 2 p Sm is determined in the same set. d n ( n is known 
because f„ is known in a neighborhood of dM. So we find that d 2 p Sm (x) is also 
determined in H m . Hence all coefficients of (|2.7j) are in fact uniquely determined 
in Ti m . 

To recover higher order normal derivatives of g and £, we continue this process, 
applying d n to equation (|2.6[) to obtain 

(2.10) (d 2 n g af3 ){v Sm ) a {v Sm )p +Hd n g a ?)d n (v Sm ) a (v Sm ) p + 

2ff Q,3 (^Ca - fl£0aPm)(«.J/3 + 1g aP d n {v Sm ) a d n {v Sm ) l3 + 

2(d n (v Sm ) n ) 2 + 2{v Sm ) n (d 2 n Cn - df lPm ) = 0. 

The second, fourth, and fifth terms in this sum have already been shown to be 
determined in H m . Therefore, so is the combined sum of the other terms: 

d 2 n g Qp (vs m Uv s Jp + 2g a0 (d 2 n ( a - d 2 n d aPm ){v Sm )p 

+ ^s m )n{dlCn-dl Pm ). 

We have already shown that d 2 p m is determined in TL m . Therefore, by differentia- 
tion, so is d 2 d a p m . By arguments similar to those of (|2.8|) and (|2.9[) . we find that 
at xq, as to — > oo, 

2(v Sm ) n (d 2 n Cn-dlp m )^Q. 

Therefore, we find that the following quantity is determined at Xq: 

d 2 n g aP (voUvo)p + (2g a Pd 2 ( a )(vo)p. 

As in the previous iteration, we regard this expression as a quadratic function in 
v. By appealing to the open nature of condition A, we find that the coefficients 
must be uniquely determined. This is true at xq, but again since condition A is 
open, it is also true in a neighborhood of Xq. Hence d 2 g a @ and d 2 8 a arc uniquely 
determined on the boundary near x . 
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Returning to equation (|2.10|) , we now find that all the remaining unknown terms 
are in fact uniquely determined in TL m - In particular, since O^Cn was already 
uniquely determined by lemma [21 d^pm must be determined in TL m . 

To determine the higher order derivatives, we successively apply d n to equation 
(|2.10p . At each stage, after discarding all terms with coefficients that are known in 
TL m , we obtain 

d k n g a0 (is m ) a (vsjp + 2g^(d^c a - d k n d aPm ){v Sm ) 

+2{v Sm ) n {d k n U~d k n +1 Pm)- 

Following (|2.8| and (|2.9|) . we find that at xo the third term goes to zero as m — > oo. 
From the last step in the previous iteration, d k p m is known in 7i m - Therefore, by 
differentiation, d k d a p m is known at xo in the limit as m — > oo. The remaining two 
terms are then regarded as a quadratic function in v and the coefficients d\g a ^ and 
d k ( a are recovered as in the previous iterations. This allows us to go back and 
determine the coefficients of the third term including d k+1 p m . Hence we find that 
the full jets of g and £ are completely determined at xq and even in some small 
neighborhood about xq. Since = g?£, the same is true for the jet of f2. 

The assumption that ui* ^ vq was used implicitly in (|2.8[) . In the case that 
w* = vq, the argument falls apart because the length of the curve approaches as 
s — ► and, consequently, d\p(x§, y Sm ) is not uniformly bounded with respect to m. 
Hence (|2.9|) does not necessarily hold and the quadratic function F(v) cannot be 
isolated. 

Case II: w* = Vo 

Lemma 4. For s > sufficiently small, p(xo,y s ) is uniquely determined by the 
scattering data. 

Proof. Recall, for small s, y s is the first intersection of j Vs with dM. Let J be a 
small interval (0, S). We shall show that if S is sufficiently small, then y s is a smooth 
curve in dM for s 6 J. This will hold for example if dM is strictly convex with 
respect to magnetic geodesies in the direction of vq. Strict convexity holds when 

(2.11) II Xo {vo)-Y(x )v 

is positive, where 77 denotes the second fundamental form of dM with respect 
to g, and Y denotes the (1, l)-tcnsor determined by the magnetic field. However, 
strict convexity is not necessary. By Taylor's formula, it suffices to show that the 
first nonzero derivative of (|2.1ip in the direction of is positive. Note that non- 
trapping and analyticity preclude the possibility that all derivatives be zero, for 
analytic continuation would then imply that the magnetic geodesic j Vo is trapped 
in dM. If the first nonzero derivative of (|2.11|) were negative, then the travel time 
necessary to reach y s from xq would be bounded from below. Hence it would be 
impossible for a sequence (y s ,w s ) to approach (xo,vq) as s — > without j Vo being 
trapped. 

We conclude that y s is a smooth curve in dM for s £ J. Therefore, p(xo,y s ) is 
a smooth function on J. By equation (|2.2j) . 

^p(zo,y s ) = {w s ,d s y s )g - (Q s ,d s y s ), 

where w s is the vector tangent to 7^ at y s . All terms on the right side are deter- 
mined by the scattering data, therefore, by integration, p(xo,y s ) is determined up 
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to a constant. By continuity at s = 0, we find that the constant must be 0, since 
lim p(xo,y s ) = 0. □ 

s— >0 

On M[ for i = 1,2, consider the quadratic function Fi(x)v = \v\g. — (Q(x),v). 
Omitting the subscript i, in boundary normal coordinates Fi takes the form: 

Fi(x)v = g af3 (x)v a vP + (v n ) 2 - ^ a {x)v a - Cn{x)v n . 

(gi, ^i) and (172, C2) have the same jet at xq if and only if 

(2.12) d^FMv = d*F 2 (x )v 

for all k > and all v in a small conic neighborhood of vo in S XQ dM . 

Suppose equation (|2.12[) fails for some t>i € S XQ dM close to «o- Then by con- 
sidering the Taylor expansion, there exists a neighborhood of (xq,v) for which 
Fi(x)v > F2(x)v for all x with a;™ > 0. Let ("f Vs )i denote the magnetic geodesic in 
M[ with initial vector v s = cos(s)^i + sin(s)9„. Then we have 



pi(x ,y s )= / Fi(x)v > / F 2 (x)v > 
'(rtv.h " / (7«Ji 

F 2 (x)v = p 2 (x ,y s ). 

(7» s )2 

The second inequality follows from the fact that (7u 3 )2 is a minimizing curve for 
the action functional defined by F 2 . This contradicts the fact that pi(xo,y s ) = 
P2{x ,y s ). □ 

Remark : The analytic assumption was only used in Case II and only to establish a 
weak convexity condition. Therefore, this theorem may be applied to nonanalytic 
systems where that condition is added to the hypothesis. 

2.5. Recovering the magnetic system in a band about dAI. The mapping 
Vo : Vi ~~ * V2 given by p.2[) is a magnetic equivalence if (<7i,f2i) and (52,^2) 
have the same coefficients when expressed in their respective boundary normal 
coordinates. By Theorem [2] and condition A on the first magnetic system, the 
respective coefficients have the same jets at least at one point in each component of 
dM. By analytic continuation the coefficients must be equal on all of dM x (— e, e). 
We conclude <p$g 2 = gi and (^^2 = 

3. Theorem [3j Recovery of the travel time data 

3.1. Jacobi fields and conjugate points. Given a magnetic geodesic 7 : [0, T] — ► 
M, let A be the operator on smooth vector fields Z along 7 defined by 

A(Z) = Z" + i?( 7 ', Z)i - Y{Z') - (V^F)( 7 ')- 

A vector field J along 7 is said to be a magnetic Jacobi field if it satisfies the 
equations 

(3.1) A(J) = 
and 

(3.2) (J',l') = 0. 
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A magnetic Jacobi field along a magnetic geodesic 7 is uniquely determined by 
specifying J and J' at a point. To see this, consider the orthonormal basis defined 
by extending an orthonormal basis ei, . . . , e n at 7(0) by requiring that 

(3.3) e' t = Y(ei) 

along 7. This extension gives an orthonormal basis at each point since 

^{e i ,e j ) = (Y(e i ),e j ) + (e i ,Y(e j )) = 0. 
Using this basis, J = ^2 fi e i an d we can write equation [ITU as the system 

n n 
i=l i=l 

where = (Y(ei),ej) and 

ay = (V r r(e ?; ) + R{i,ei)i - Y(Y(e l )) - (V ei F)( 7 '), ej) . 

This is a linear second order system, and therefore it has a unique solution for 
each set of initial conditions. Moreover, since the metric and magnetic fields are 
analytic, so will be the solutions. 

It is easy to see from the definition that 7' is always a magnetic Jacobi field. 
Unlike the case of straight geodesies, this is the only magnetic Jacobi field parallel to 
7'. Another difference from the straight geodesic case is that magnetic Jacobi fields 
that are perpendicular to 7' at t = do not stay perpendicular for all t. For this 
reason we will sometimes consider instead the orthogonal projection J = J — f-f' 
where / = (J, 7'}. The parallel component f^ 1 is uniquely determined by and 
J(0) since 

/' = (J', 7') + (J, 7") = (J,Y( 7 ')) = (J X ,Y( 7 ')) . 

Magnetic Jacobi fields are the variational fields corresponding to variations through 
magnetic geodesies. This can be seen by considering the variation 

f(t,s) = ls (t) = ex P ^ s) (t6(s)) 

where r(s) is any curve with r'(0) = J(0) and 6(s) is a vector field along t with 
9(0) = 7'(0) and 0'(0) = J'(0). It is a straight forward computation to check that 
the variational field 0) is a magnetic Jacobi field (see [He]). 

There is a close relation between magnetic Jacobi fields and conjugate points, 
analogous to the straight geodesic case. This can be summarized in the following 
Proposition, proved in [He] . 

Proposition 3.1. Let j$ : [0, T] — > M be the magnetic geodesic with 7(0) = x. 
The point p = 7(^0) is conjugate to x along 7 if and only if there exist a magnetic 
Jacobi field J along 7, not identically zero, with J(0) = and J (to) parallel to 7'. 

Moreover, the multiplicity of p as a conjugate point is equal to the number of 
linearly independent such Jacobi fields. 
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3.2. Proof of Theorem [3l Consider magnetic geodesies 71 : [0, Ti] — > Mi and 

72 : [0, T2] — > Af2 that satisfy condition B and such that 72(f) = ifio (71(f)) for small 
f. Suppose, moreover, that 71 is never tangent to dM±. Since both systems have 
the same scattering data they must exit the manifold at corresponding points and 
directions, that is 72(^2) = Vo(7i(Ii)) and 72CT2) = fo*(l'i(Ti)). We will call such 
geodesies correspondent, and, through ipa, shall implicitly identify 71 (0) to 72(0) as 
well as their neighborhoods and tangent spaces. 

We want to compare variations of 71 and 72 by defining for a variation 

f 1 (t,s) = exp* (s) (te(s)) 

of 71 a corresponding variation 

f 2 (t,s) = exp^ o(T{s)) {t<p {6{s))) 

of 72. Nonetheless, fi(T\,s) will be in V for small s and that is enough for our 
purpose. Therefore /i(Ti, s) = / 2 (T(s), s) where T(s) : (— e, e) — > R is diffcrentiable 
with T(0) = T 2 . By taking derivatives with respect to s we see that 

^m,0) = ^(T2,0)+T'(0)72(T 2 ), 
os os 

so corresponding Jacobi fields agree at their endpoints up to a term tangent to 7-. 

We want to consider all Jacobi fields along 71 with J(2i) = 0. For this we 
define the magnetic Jacobi tensor J\ with the initial conditions Ji{T\) = and 
J[(Ti) = Id. This is a (1, 1) tensor in the space of vectors orthogonal to r y[. For 
any vector 77 at 71 (Ti ) perpendicular to j[ , the vector field J (t ) = J\ (77) (t) is the 
perpendicular part of the Jacobi field with J(0) = and J'(0) = ?/. 

Correspondingly, define the magnetic Jacobi tensor J 2 with Ji(T-i) = and 
•^2(^2) = Since correspondent Jacobi fields agree at their endpoints, up to a 
parallel component, so will the corresponding Jacobi tensors. Therefore 3\{t) = 
J~2(t) for small t. Moreover, this is an analytic function and we can use analytic 
continuation to see that they must agree for all t , as long as both are well defined. 

Assume, without loss of generality, that T\ < Ti. Then, Ji{i) = J<z(t) for 
< t < Ti. In particular, at T\ we have JiiTi) = Ji{T\) = 0. which would imply 
that 72 (Ti) is conjugate to 72^2), of order n — 1. Since 71 satisfies condition B 
this can't be the case, therefore T2 = Ti. 

To see that the length of all magnetic geodesies agree, note that magnetic 
geodesies that are tangent to dM form the boundary of Q. In particular they 
are limit points, and form a set of measure 0. Therefore, we know that the lengths 
of 71 and 72 agree for a dense subset of Q\. The length of geodesies is a continuous 
function on Q, so the difference in lengths is also a continuous function that is in 
a dense subset. This implies that the lengths agree for all magnetic geodesies. 

3.3. About the necessity of condition B. In the presence of conjugate points 
of order n — 1 the proof above would say that the lengths may differ by the dis- 
tance between conjugate points of such order, and that Ji(T\ — s) = J\{T\ — s) = 
tp Ji^Ti — s), where ip is a local isomctry of TM2 allowing the possibility that the 
diffeomorphism tpo does not extend in a compatible way through 71 and 72 . So Ji 
has a period, up to isometry, of length T2 — Ti. 
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In the case of straight geodesies (f2 = 0) we can see this behavior by looking at 
the following examples. 

Let Mi be an analytic subset of the sphere S 2 that contains the southern hemi- 
sphere. Let Mi be a region of WP 2 with the same boundary. We can see it as 
removing the southern hemisphere from M\ and identifying antipodal points in the 
equator (See figure [l}. This manifolds have the same scattering data but any geo- 
desic in Mi that passes through the southern hemisphere will be longer by it than 
the corresponding geodesic in M 2 . 

Another similar example is when Ni is a subset of the sphere S 2 that contains 
both poles and a meridian r. To build N2 we cut open along r and glue another 
S 2 opened at the meridian in such a way that whenever you cross r you move from 
one sphere to the other (see figure [2]) . These manifolds have the same scattering 
data but any geodesic in N\ that passes through r will be shorter by 2tt than the 
corresponding geodesic in iVjj. This is not a Ricmannian manifold, since the poles 
are singular points. But it is orientable and the Jacobi tensors will be periodic with 
a period of 2ir. 

The need of condition B is clear in the examples above, where we can see that 
the presence of a conjugate point of order n — 1 allows the length of a geodesic to 
change. The question remains if it is necessary as an independent hypothesis. It 
seems to be that, like in the examples above, conjugate points of order n — 1 only 
appear when there are trapped geodesies. If this is true, it would make condition 
B empty in this setting. 

4. Theorem [4] Construction of the magnetic equivalence 

We have two non-trapping, analytic magnetic systems (Mi,gi, f2j) with the same 
scattering data and travel time data. Moreover, the mapping ipg : V\ — > Vi is 
a magnetic equivalence. To prove Theorem 2J we must extend (po to a magnetic 
equivalence <p : M[ — > M' 2 . 

This has been done in the absence of magnetic fields in [Vj . The proof carries 
over in the presence of magnetic fields with only changes in notation. Therefore, 
we shall only sketch the proof here. 




Mi 



M 2 



Figure 1. 



4.1. Extension of <pg to M{. 
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Figure 2. 

Proof. First we define a mapping dp : SMi — ► M2. We then show that all vectors 
lying over the same base point in Mi map to the same point in M2, and that the 
induced map tp : M\ — > M2 agrees with ipo on their common domain. 

Let (x,6) G 5Mi, and consider the magnetic geodesic "f x ,e{t)- Going backwards 
along 7a;, e, we will eventually leave M\ by the non-trapping assumption. We define 

T = T (x,6) = mf{t > : j(x,9)(-t) t Ah}. 

Tq > satisfies ^f Xt e(— Tq) G 9A/. Moreover, by the fact that dM is analytic, it is 
not hard to see that j Xt g (—To — s) £ V\ \ Mi for all s such that < s < S for some 
sufficiently small S = 5{x,9). We arbitrarily choose 

-T = -T(x,6) g (-T -<5,-T ), 
and let (^a;^, G SM[ be given by 

(zx,e,£x,e) = {lx,e{-T),j' xe (-T)). 
Let (y x ,e,rix,e) G be its image by tp , 

{yx,d,rjx,d) = <po*{zx,e,tix,e)- 

Define 

(p{x,6) = 7fe j9 ,^,9( T )- 
Proposition 2. (p is a well defined function on SM\ with values in Mi- 
Proof. It must be shown that the value of <p(x, 6) is independent of the choice of 
T G (—T — (5, —T ), and that the values do belong to M 2 . The first statement 
follows from the fact that the two magnetic systems are equivalent via ipo in V. 
The second follows from the fact that the two systems have share the same length 
data. See [V] for details. □ 

Proposition 3. The value of<p(x,6) is independent of the direction 9. 

Lemma 5. Let M' be an open manifold with analytic metric g. Then for every 
xq G M' , there exists a positive number r such that the squared distance function 
is analytic on the set 

A r (xo) = {(x,y) : d(x,x ) < r,d(x,y) < r}. 
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If K is a compact set contained within the interior of M' , then there is an open 
O C M' containing K and a positive number r such that the squared distance 
function is analytic on the set 

AoA K ) = {(x,y) ■ x £ 0,d(x,y) < r}. 

See jV] for proof. 

Proof. Since the sphere \8\ g = 1 is connected, it is sufficient to show that for fixed 
x, <p(x, 9) is locally constant in 9. 

Fix (xo,#o)j an d let N C S Xo Mi be a neighborhood of 9$. Let (zo^Co) and 
(UOtVo) correspond to (xq,9q) following the notation above. Then we have 

x =7* ,£o( T ) 

<p(x ,e ) =7 VO ,r,o(T). 

For all 9 £ N, let (z,$) = (l xo A-T),j' Xo , g (-T)), so that 

xq = lz,d T )- 

Let (y, 7]) = (fio*(z, £). It can be shown that if TV is sufficiently small, 
(4.1) <p(x ,6)=~/ Vtri {T). 

Equation (|4.ip seems to be a direct application of the definition of (p but in fact, 
between xq and z, the curve 7k 0i 6( can pass across the boundary in and out of M\. 
Therefore, T is not necessarily a valid choice for T(xq, 8) according to the definition. 
Nevertheless, using the fact that the two magnetic systems have the same scattering 
and length data, equation (|4.1|) still holds. 
According to lemma the functions 

Pi(t) =^ 1 (7«o,e (*) 5 7*,c(*)) ) and 

will be analytic wherever their values are less than some fixed constant r. 

Since ipo is a magnetic equivalence, it carries magnetic geodesies into magnetic 
geodesies. (y,rj) = ipo*{z^). Therefore, for t near 0, we have 

Po(7*,£(*)) = Iv.vi*)- 

A magnetic equivalence is, in particular, an isomctry of metrics, so we conclude 
that for t near 0, pi(t) = P2(t). By analytic continuation, the two function must 
be equal for all t up to T. Therefore, P2(T) = pi(T) = 0. This implies that 

<p(x ,6) = j y ,r,( T ) = ly ,vo( T ) = <p(x o ,0 o ). 

□ 

We define (p : M{ -> M' 2 by 

x^(p(x,ff)i f° r x £ Mi 

x *-^(fio (x) , for x £ M[ \ Mi . 

The magnetic geodesic flows on the respective sphere bundles of M\ and M 2 are an- 
alytic. This implies that ip is an analytic mapping. The fact that flvx = follows 
directly from the definition of tp and the fact that tpo is a magnetic equivalence. By 
analytic continuation, tp must satisfy 

V>*92=9x, (p*Q.2 = U\. 
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